Non-Markovian qubit decoherence during dispersive readout 
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We study qubit decoherence under generalized dispersive readout, i.e., we investigate a qubit 
coupled to a resonantly driven dissipative harmonic oscillator. We provide a complete picture 
by allowing for arbitrarily large qubit-oscillator detuning and by considering also a coupling to 
the square of the oscillator coordinate, which is relevant for flux qubits. Analytical results for 
the decoherence time are obtained by a transformation of the qubit-oscillator Hamiltonian to the 
dispersive frame and a subsequent master equation treatment beyond the Markov limit. We predict 
a crossover from Markovian decay to a decay with Gaussian shape. Our results are corroborated by 
the numerical solution of the full qubit-oscillator master equation in the original frame. 

PACS numbers: 03.65.Yz, 03.67.Lx, 42.50.Dv, 85.25.Cp 
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I. INTRODUCTION 

The final readout of the qubit statepresents an essen- 
tial part of any quantum algorithm [IJ. For solid state 
qubits, it is typically realized by coupling the qubit to 
a harmonic oscillator such that the oscillator frequency 
undergoes a shift whose sign depends on the state of the 
qubit. This shift can be probed by driving the oscil- 
lator at its bare frequency, with the consequence that 
the phase of the response provides information about the 
qubit state 0-0] • Typically, one works with a qubit- 
oscillator detuning that is slightly larger than the respec- 
tive mutual couplingwhile still being much smaller than 
the qubit splitting 0, H[ . Then the resulting frequency 
shift of the oscillator can be derived from the Rabi Hamil- 
tonian via a transformation to the so-called dispersive 
frame Q. For even stronger detuning, the transforma- 
tion becomes more involved, but nevertheless the sign of 
the dispersive shift depends on the qubit state @ and, 
thus, qubit readout remains possible. For an extreme 
detuning, such that the oscillator frequency exceeds the 
qubit splitting by far, a measurement protocol has been 
proposed by which one can reconstruct information about 
coherent qubit oscillations from recorded data [10]. Even 
though similar readout is possible by driving the qubit 
directly [III EH , the oscillator plays a constructive role 
as band pass. 

For driving the oscillator and for measuring its re- 
sponse, the setup must be coupled to auxiliary electronic 
circuitry, which represents an environment that eventu- 
ally destroys the phase of the qubit. Generally informa- 
tion about the qubit state can be obtained only at the 
rate at which the qubit coherence decays [l3j]. Thus a 
quantitative understanding of qubit decoherence stem- 
ming from the coupling to a resonantly driven oscillator 
is inevitable for the design of dispersive readout schemes. 
For sufficiently small detuning, such that the rotating- 
wave approximation underlying the Rabi Hamiltonian 
holds, the decoherence rate follows from an intuitive con- 
sideration in which the shot noise of the cavity photons 
randomizes the qubit phase Q. This result will emerge 



as limiting case of our more general picture. 

For weakly dissipative quantum systems, Bloch- 
Redfield theory 0, [l5| represents a natural framework 
for studying decoherence, in particular when memory ef- 
fects are minor. When qubit decoherence stems from 
a dissipative harmonic oscillator, however, its naive ap- 
plication may significantly overestimate the decoherence 
rate [H, [TtJ because peaks in the effective spectral den- 
sity of the oscillator |18l - l20l ] cause non-Markovian behav- 
ior. This may in particular be the case for the ultra- 
strong qubit-oscillator coupling which marks a recent 
trend ^MM- 

Here we present a global picture of the qubit deco- 
herence during dispersive readout. We consider both 
linear and quadratic qubit-oscillator coupling as well as 
arbitrarily large detuning, while we do not account for 
higher-order corrections to the dispersive shift 25| and 
non-linearities in the oscillator potential (26|. While 
in Ref. [l9[ this problem was studied for an oscillator 
at thermal equilibrium, we focus on the semiclassical 
regime, in which the oscillator is governed by the driv- 
ing. In Sec. [Til we introduce our quantum master equa- 
tion for the full qubit-oscillator-bath model. Section IIIII 
is devoted to an analytical derivation of the dephasing 
time. We develop a picture in which the driven dissi- 
pative oscillator acts as a bath which is eliminated in 
second-order perturbation theory but beyond a Markov 
approximation. In order to corroborate the resulting de- 
phasing time, we numerically solve in Sec. HVI the quan- 
tum master equation for the qubit and the oscillator in 
the original frame. 



II. QUBIT-OSCILLATOR MODEL 

The qubit coupled to the oscillator is described by the 
Hamiltonian 

Hq = —a z +hgia x (a^ +a)+hg 2 cr x (a^ +a) 2 +hujQa ! a, (1) 
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where g\ is the strength of the dipole interaction which 



is linear in the oscillator displacement a 



The cou- 



pling to the square of the displacement with strength g 2 
is mainly relevant for flux qubits and can be controlled 
to some extent [271.128}. 

The system state can be probed via a coupling an 
external circuitry which we model by the system-bath 
Hamiltonian 



fli 



bat h 



- (ot + cOVct&t + b v ) + yinu v blb v , (2) 



where b v is the annihilation operator of a bath mode v. 
The influence of the bath is determined by the spectral 
density I(oj) = 27r c£<J(c<; — w v ), which we assume to 
be ohmic, i.e., I(oj) = "/lo/ujq with the oscillator damping 
rate 7. 

An external ac driving corresponds to one particular 
bath mode being in a highly excited coherent state. It 
may be described as classical oscillation, such that the 
oscillator experiences a driving force Acos(ttt — <po). This 
corresponds to the driving Hamiltonian 



H F (t) = Axcos(Vtt 



(3) 



For convenience, we have introduced the dimensionless 
position and momentum operators x — {a) + a)/\/2 and 
p = i(at — a)/v2j respectively. 

Within the usual Born-Markov approximation for the 
oscillator (see Appendix^}, one obtains the master equa- 
tion 



R = - \ [Ho + H F (t)\ - ^[x, [p, R]+] 



(4) 



for the joint density operator R of the qubit and the oscil- 
lator. It provides all numerical results presented below. 
In case of low temperatures, fc^T <C ?u>Joi the hyperbolic 
cotangent assumes a value close to unity. Notice that in 
contrast to a parametric driving, the linearly coupled ac 
force does not affect the dissipative terms [29| . Moreover, 
we assume that the qubit couples only weakly to the os- 
cillator, so that its influence on the dissipative kernel of 
the oscillator is negligible. 



III. ANALYTICAL ESTIMATE OF THE QUBIT 
DEPHASING TIME 



Extracting from the master equation (j4]) an analyti- 
cal expression for the qubit decoherence time represents 
a formidable task. Thus, we have to rely on several 
approximations that make use of the conditions under 
which dispersive readout may be performed. We start by 
a transformation to the dispersive frame which yields a 
coupling to the square of the position coordinate of the 
driven oscillator. The relevant influence on the qubit is 



determined by the corresponding auto correlation func- 
tion which we evaluate in the semiclassical limit. It be- 
comes stationary only after averaging within a rotating- 
wave approximation over the initial phase (f>o of the driv- 
ing. Finally, the resulting non-Markovian master equa- 
tion for the qubit is solved for short times. 



A. Transformation to the dispersive frame 

For the discussion of dispersive readout, the disper- 
sive picture of the qubit-oscillator Hamiltonian Hq has 
proven useful @, [l(| HH . In order to capture also very 
large qubit-oscillator detuning, we need to perform the 
according transformation beyond rotating-wave approxi- 
mation. This yields the effective qubit-oscillator Hamil- 
tonian M 



Hq = yCT z + fi(A||(T a + \±<7 X )X 2 + ^-{x 2 + p 2 ) (5) 



with the coupling constants 



\ ffi ffi 

Am = 1 . 

e — u!q e + ojQ 
X± = 2.92- 



(6) 
(7) 



Thus, the coupling linear in the oscillator coordinate has 
turned into a quadratic coupling with strength Aii, while 
Aj_ has been introduced for a unified notation. In cor- 
respondence to the orientation of the coupling operators 
on the Bloch sphere, we refer to the coupling terms as 
"longitudinal" and "transverse" , respectively. 

If only energy-conserving terms in the qubit-oscillator 
coupling were considered [8j, Aii would be given by only 
the first term of Eq. ([6]) and, thus, be cx (e — wo) -1 . By 
contrast, the second term of Eq. ^ turns the frequency 



dependence into Ay 



This means that for 



positive detuning, ujq > e, the counter-rotating terms di- 
minish the dispersive shift. Since we will find that deco- 
herence grows with A|| , the coherence time is larger than 
predicted within rotating-wave approximation 8] . 

The interpretation of the effective interaction is that it 
shifts the oscillator frequency by ±(A| + A^) 1 / 2 , where 
the sign depends on the state of the qubit. There- 
fore, probing the oscillator frequency provides informa- 
tion about the latter. In this work we are interested 
in the qubit decoherence that stems from this coupling. 
When transforming to the dispersive frame, the dissipa- 
tive terms of the master equation ^ acquire a contri- 
bution from qubit operators [30l. [3l|. which, however, is 
negligible as compared to the terms considered below. 



B. Driven oscillator as effective bath 

We now treat the oscillator as environment coupled to 
the qubit coordinate Y = A||<7 2 + Xj_a x via the Hamil- 
tonian H lnt — frYrj with the environmental fluctuations 
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r\ determined by the operator x 2 . Its expectation value 
(x 2 ) yields a correction to the qubit Hamiltonian of the 
order A|| . Therefore, the impact of (x 2 ) on the dissipative 
terms is already beyond the order considered herein and 
can be omitted, such that the relevant fluctuation reads 

i 1 = x 2 -{x 2 ). (8) 

We assume weak dissipation such that the bath can be 
eliminated within second order perturbation theory. This 
is in accordance with our observation of predominantly 
coherent time evolution; see the numerical results in 
Sec. IIVI Then the dissipative part of the master equa- 
tion for the qubit density operator p reads 

p = - [ dt{s vv (t,t')[Y,[Y(t-t'),p}} 

Jo 1 (9) 

+iA vr ,(t,t')[Y,[Y(t-t'),p] + }}. 

It is non-Markovian owing to its explicit time-dependence 
and the corresponding lack of a semi-group property 
Here [ , ] + denotes the anti-commutator, while 

Y(t) = \\a z + \±{a x cos(ei) - a y sin(rf)} (10) 

is the qubit part of the coupling in the interaction picture. 
3^ and A m are the real part and the imaginary part, 
respectively, of the effective bath correlation function 

C nn (t,t') = (v(tW)) = S, m (t,t')+iA vv {t,t'), (11) 

which is not time homogeneous due to the driving. The 
evaluation of S vri in the limit in which dispersive readout 
is performed is a cornerstone of our analytical treatment. 

Since we consider measurement schemes that rely on 
the response to deterministic driving, the fluctuations are 
small so that we can linearize in the oscillator position 
fluctuation 8x = x — x{t) and work with the approxima- 
tion 

r}{t) = 2x(t)5x(t). (12) 

Then the auto correlation function of the effective bath 
coordinate r\ becomes 

C m {t,t') = Ax{t)x{t') (5x(t)Sx(t')}, (13) 

where the term with the angular brackets is the position- 
position correlation function Cg x g x of the dissipative har- 
monic oscillator. Owing to the linearity of the oscillator's 
equation of motion, it is independent of the driving and, 
thus, stationary. 

The response to the classical driving can be ex- 
pressed in terms of the oscillator Green's function G(w), 
Eq. ([£3l> . This yields 

x(t) = VMcos(Qt - 0o - 4>), (14) 

where h = iyl 2 |G(c<j)| 2 is the mean cavity photon num- 
ber and <f> is the phase of the Green's function, while 



4>q is the unknown initial phase of the driving. For a 
harmonic oscillator in the weak coupling regime, the re- 
sponse to the fluctuations of the external circuitry is con- 
veniently com put ed with the help of the quantum regres- 
sion theorem [32j , see Appendix IA 31 For low tempera- 
tures, ksT <C fiwoi it can be approximated by 

Ssx8x(t) = ie-^/ 2 cos(w r), (15) 

where r = t — t' . Equations (|Ti|) and (fT5")) allow us to 
evaluate the correlation function (fT3")) . After performing 
an average over the initial phase <pQ, we obtain for its 
symmetric part the time homogeneous expression 

S vv {t) = 27le~ 7T/2 cos(f7r)cos(w T). (16) 

The phase average represents a rotating-wave approxi- 
mation and is possible since qubit decoherence and dis- 
sipation are much slower than the coherent oscillator dy- 
namics. The correlation function flT6"|) possesses four res- 
onance peaks of width 7 at the frequencies u = ±ujo ± f2, 
where for resonant driving, Q = u>q, the two central peaks 
coincide at zero frequency. 

The peaks of the spectral density correspond to long- 
time correlations of the quantum noise that may lead 
to non-Markov effects. Therefore, a treatment with a 
fully Markovian master equation is not appropriate [l7| . 
We thus employ ideas that have been used for studying 
decoherence due to 1// noise [33T4371 ] . We analyze the 
longitudinal and the transverse dephasing separately by 
setting either g\ or to zero. 

C. Coherence decay for linear qubit-oscillator 
coupling 

We assume that the qubit is initially in the state 
(It) + |4-))/v2) i-e., in a coherent superposition of the 
eigenstates of the qubit Hamiltonian (He/2)a z . Then the 
corresponding off-diagonal elements of the density matrix 
in the eigenbasis are both 1 /2 and undergo an oscillatory 
decay, = |exp{— iet — A(t)}, where in the Markov 
limit, A(t) = Tt. It is straightforward to demonstrate 
that then the purity P = tr p 2 , being our measure of 
coherence, evolves as 

P{t)= l -[l + e- 2 ^)]^l~K{t). (17) 

The approximation holds for short times at which the 
purity still lies significantly above 1/2. A typical time 
evolution is depicted in Fig. [1] below. It demonstrates 
that the purity decay is indeed not necessarily a simple 
exponential, but may have Gaussian shape. 

The still unknown function A(t) will be determined 
from the master equation ([9]) for the density matrix ele- 
ment p^i at short times yielding for 52 = 

A|,(t) = -£± = r|,(t) (is) 
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with the time-dependent decoherence rate 



r,|(f)=4A? / drS m {T). 



(19) 



The index || refers to the longitudinal qubit-oscillator 
coupling in the dispersive Hamiltonian ([5]). Notice that 
in the original Hamiltonian (JlJ, this coupling is trans- 
verse. In the following, we evaluate this rate for resonant 
driving and weak oscillator damping, 7 -C ujq = 17. 

Inserting the effective spectral density (fT6|) into 
Eq. (H11) yields 



8nA ; 



7 



(20) 



where we have neglected terms oscillating rapidly with 
frequency 2wo- By straightforward time integration, we 
obtain 



8nA 2 



A ll( t ) = -^(^ + 2e- 7t/2 -2) 
'2nA|i 2 for7t<l, 
in\f,t/j for jt » 1. 



(21) 
(22) 



Inserting this approximation into our ansatz for p^i re- 
veals that during an initial stage, the coherence decays 
like a Gaussian p^ ~ exp(— t 2 /Tq) with the time scale 



T G = 



2n\l 



(23) 



Thereafter, normal exponential decay p^i ~ exp(— t/Tj\/) 
sets in, where 



7 



8nA 2 



(24) 



Since both approximations in Eq. (|22p are never smaller 
than the exact expression, we connect the two limits by 
choosing at each time the smaller one, i.e., 



Ay(i) ^mm(t 2 /T 2 ,t/T M ). 



(25) 



This implies a crossover from Gaussian to Markovian de- 
cay at time t c = 4/7. 

As criterion for "significant dephasing", we use A > 
1/4, which means that the off-diagonal matrix element 
Pfl has decayed by at least 22 percent. For larger val- 
ues of A, the visibility of coherent oscillations is already 
quite small. Therefore, the relevant dephasing is Gaus- 
sian if 1%/Tq > 1/4 or, equivalently, 128nA| > 7 2 . In the 
opposite case, the Gaussian stage can be ignored and co- 
herence fades away during time Tm/4. In combination, 
this yields the dephasing time 



W4=^Af forl28rlA|< 7 2 , 
T' ! , „ (26) 



T G /2 = -4— for 128nA 2 > 7 2 . 

V 8fiA ii 11 



The first line holds for the Markovian behavior found for 
weak coupling. Notice that A|| is an effective coupling 
constant that becomes smaller with increasing detuning 
I e — ojq j . Thus, for small detuning (but still within the 
dispersive limit) and for large photon number, we expect 
Gaussian decay. 

At this stage, it is interesting to establish a connection 
to Refs. [8|,|l3(, where the fluctuation of the cavity photon 
number leads to a fluctuating qubit splitting and, thus, 
randomizes the qubit phase. Then, for |A||| < 7 and 
A^ = 0, one finds that the off-diagonal matrix elements 
of the density operator decay with a rate = 1/Tm — 
8A 2 n/ 7 HG1, which 

is accordance with our result in the 

Markov limit. 



D. Coherence decay for quadratic qubit-oscillator 
coupling 

For <?2 7^ 0, the situation becomes considerably more 
complicated, because the dissipative terms in the master 
equation @ couple all density matrix elements to each 
other. Therefore, one can no longer obtain a closed first- 
order equation for p^ such as Eq. (fl"8j) . For this reason 
we attempt an analytical solution only in the Marko- 
vian regime. In doing so, we perform the time inte- 
gral in Eq. Q until infinity such that we obtain a time- 
independent Bloch-Redfield master equation. The deco- 
herence rate is conveniently extracted from the equiva- 
lent equation of motion for the Bloch vector s = (a). By 
straightforward algebra the latter emerges as ds/dt = 
Ms + h with the dynamical matrix 



M 



and the decay rate 




A9 5 



2 U VV 



7 2 + 4e 2 ' 



(27) 



(28) 



The inhomogeneity h stems from the second line of 
Eq. (j9|) and determines the stationary state which is not 
relevant in the present context. 

We proceed by computing the eigenvalues of M to low- 
est order in the dissipative terms, which yields — 2r^ and 
±ie — r_L (38|. The latter correspond to the decaying os- 
cillations of p^i, which reveals that the transverse deco- 
herence is determined by Aj_ = T±t. Thus, we obtain in 
the Markov limit the time scale T^j = 1/T± and, thus, 
the dephasing time 



rr 



1 



_ 7 2 + 4e 2 

4r ± ~~ 167nA 2 L 



(29) 



This result holds under two conditions. First, it is re- 
quired that the time-integration in the master equa- 
tion §§§ can be extended to infinity, which is possi- 
ble if the decay time of the effective bath correlation 
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function (|T6]) is much shorter than the dephasing time, 
I/7 < T* , which means 16nA^ < 7 2 + 4e 2 . Second, the 
qubit frequency must be within the oscillator linewidth, 
i.e., e < 7, because otherwise the oscillator would shield 
the qubit from the external circuitry. Then higher-order 
processes may dominate while the master equation (|9|) is 
of only second order. The latter condition is also essential 
for an application that we have in mind, namely time- 
dependent dispersive qubit readout via a high-frequency 
oscillator [uj. 



IV. NUMERICAL DETERMINATION OF THE 
QUBIT DEPHASING TIME 

The numerical computation of the dephasing time T* 
from the full master equation Q is possible only in a 
restricted parameter regime for various reasons. First, 
at resonant driving, the stationary state of the oscillator 
has mean photon number n = A 2 /2-f 2 . The assump- 
tion of the oscillator being in its semiclassical limit is 
fulfilled only for n 3> 1 or equivalently A 3> 7. Second, 
our considerations in Sec. [TT] require that the oscillator 
reaches its stationary state during a time much shorter 
than the qubit dephasing time, i.e., for 1/7 <C Tj|,Tj_. 
Finally, the stationary photon number is limited by com- 
putation time, which grows with the number of oscillator 
Fock states needed for numerical convergence. While we 
find a good agreement of the numerical and the analyti- 
cal results already for n > 5, the natural expectation is 
that the agreement even increases with the mean pho- 
ton number, because then the oscillator becomes more 
classical. 

We start with the oscillator in the coherent state that 
corresponds to the stationary classical solution in the ab- 
sence of the qubit. For the qubit itself, we use as initial 
state the superposition (||) + \\))/^/2. The dissipative 
time evolution of both the qubit and the oscillator is ob- 
tained by numerical integration of the full master equa- 
tion with the original Hamiltonian ([1]). Thus we im- 
plicitly also test the validity of the dispersive picture in 
the presence of a heat bath. 

The central quantity of our numerical study is the 
time evolution of the purity Pit) — trp 2 (t) from which 
we determine the dephasing time T* by the criterion 
P(T*) = ^(l + e" 1 / 2 ), i.e., A(T*) = 1/4 as above. More- 
over, we use P(t) to decide whether decoherence decays 
like a simple exponential or like a Gaussian. A formal 
procedure for the distinction is fitting P(t) for short times 
to the ansatz P(t) — 1 — a^/t — a^t 2 . The decay is mainly 
Markovian or mainly Gaussian depending on which rate 
ajv/ or clq is larger and thus dominates. 



A. Linear qubit-oscillator coupling g\ 

Figure Q] depicts the time evolution of the qubit ex- 
pectation value (<r x ) which exhibits decaying oscillations 
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FIG. 1. (Color online) Typical time evolution of the qubit op- 
erator (a x ) (solid line) and the corresponding purity (dashed) 
for fi = Uq = 0.8e, gi = 0.02e, 7 = 0.02e, and driving am- 
plitude A = 0.06e such that the stationary photon number 
is TT, = 4.5. Inset: Purity decay shown in the main panel 
(dashed) compared to the decay given by Eq. (|17[) together 
with Eq. J2TJ (solid line) . 
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FIG. 2. (Color online) Dephasing time for purely linear qubit- 
oscillator coupling ((72 = 0), resonant driving, fi = cjq, and 
oscillator damping 7 = 0.02e. The amplitude A = 0.07e 
corresponds to the mean photon number n — 6.125. Filled 
symbols and dashed lines refer to predominantly Markovian 
decay, while for Gaussian decay, stroked symbols and solid 
lines are used. 



with frequency e. The parameters correspond to an inter- 
mediate regime between the Gaussian and the Markovian 
dynamics, as is visible in the inset. 

In Fig. [21 we compare the decay time Tjf with the 
analytical result (|26l) for various values of the oscillator 
damping and the qubit-oscillator coupling g\ as function 
of the oscillator frequency. Whether Markovian or Gaus- 
sian decay dominates is visualized by filled and stroked 
symbols, respectively. We have skipped the regime very 
close to resonance, |Aii| < 5|wo — e|, since there the disper- 
sive Hamiltonian ([S]) is not valid and so far no dispersive 
readout protocol has been proposed. 

Our prediction for regimes with non-Markovian decay 
is confirmed by the numerical solution rather well. No- 
tice however, that in the crossover regime, our formal 
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criterion for Markovian decay provides a unique answer, 
even though the respective other decay may already con- 
tribute significantly. The agreement of the numerically 
found border with our prediction corroborates as well 
the crossover time t c conjectured above. Concerning the 
values of Tjf , we observe a good overall agreement with 
the tendency that the analytical result slightly underesti- 
mates T|*. In the regime of Gaussian decay, a Markov ap- 
proximation would yield a significantly smaller coherence 
time. This means that, interestingly enough, the qubit 
stays coherent for a longer time than is expected from 
Bloch-Redfield theory. For large oscillator frequency, 
uj 3> e, also the predicted behavior T|* oc A^ 2 oc |e 2 — w 2 | 2 
is confirmed. This substantiates the relevance of the 
counter-rotating correction to the dispersive shift [9| . 

B. Quadratic qubit-oscillator coupling g2 

We proceed as above, but consider the coupling to the 
square of the oscillator coordinate in the Hamiltonian (fTJ) , 
while setting gi = 0. Even though the linear coupling gi 
can be controlled to some extent [13, [H| , it is probably 
hard to turn it off completely. Still our choice has rele- 
vance in the limit of large detuning in which the effective 
dispersive coupling An becomes rather small, see Eq. ([6]). 
Then for realistic values of g^ for flux qubits, a protocol 
for recording coherent time evolution has been proposed 
10]. A necessary condition for this is an oscillator band- 
width of the order of the qubit splitting, such that the 
oscillator does not filter out the information about the 
coherent qubit dynamics. Therefore, we will choose an 
oscillator with the rather large frequency Q = ujq = 5e 
and with damping up to 7 = e. 

Figure [3] shows the numerically obtained coherence 
times and whether the decay is predominantly Gaussian 
or Markovian. For the large oscillator damping 7 = 6, 
the conditions for the validity of the (Markovian) Bloch- 
Redfield equation stated at the end of Sec. IIIIDI hold. 
Then we observe a good agreement of the numerically 
obtained T£ an< ^ Eq. (|29l) . There seems a slight system- 
atic deviation for large values of 92- In this limit, how- 
ever, the numerical results are not very precise, because 
the purity falls off already during the first few oscillation 
periods. This rapid decay hinders a precise numerical de- 
termination of the dephasing rate. Thus, the agreement 
is still within the numerical precision. 

For small oscillator damping, 7 <C e, Bloch-Redfield 
theory is not applicable, as discussed above. Moreover, 
the short time evolution of the purity is already too 
complex for a comprehensive prediction of the decoher- 
ence time. Nevertheless our numerical solution provides 
some hint on the decoherence process. We find that the 
crossover from simple exponential decay to decay with 
a Gaussian shape occurs at smaller values of the qubit- 
oscillator coupling. For small (72, the qubit stays coher- 
ent slightly longer, while for large g2, coherence decays 
a bit faster as compared to the case 7 = e. In both 
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FIG. 3. (Color online) Dephasing time for purely quadratic 
qubit-oscillator coupling (g\ — 0), resonant driving at large 
frequency, Q — luq — 5e, and various values of the oscillator 
damping 7. The driving amplitude is A = 3.57, such that 
always n = 6.125. Filled symbols mark Markovian decay, 
while stroked symbols refer to Gaussian shape. The solid line 
depicts the value obtained for 7 = e in the Markov limit. The 
corresponding numerical values are connected by a dashed 
line which serves as guide to the eye. 

regimes, the 7 dependence of T±_ is weak. This disproofs 
the Markovian theory for 7 C e, because the latter pre- 
dicts Tm oc I/7, which stays in contrast to our numerical 
result. 



V. DISCUSSION AND CONCLUSIONS 

A qubit undergoing dispersive readout, i.e., one that 
is coupled to a resonantly driven dissipative harmonic 
oscillator, experiences decoherence from a rather exotic 
effective environment. The latter's main properties stem 
from the small oscillator linewidth, the strong driving, 
and a coupling coordinate that does not commute with 
the qubit Hamiltonian. Nevertheless, it has been possible 
to analytically obtain essential and concise information 
about the decoherence process. Our approach is based 
on a transformation to the dispersive frame, which turns 
the linear coupling into phase noise. In doing so it is 
crucial to perform the transformation beyond rotating- 
wave approximation, in particular for studying the far 
detuned case where the non-rotating corrections are of 
the same order as the rotating- wave terms. For the sub- 
sequent analytical treatment, we have derived the de- 
phasing time within our picture of an effective spectral 
density provided by the driven harmonic oscillator in the 
semiclassical limit. At the same time it has turned out 
that a peaked spectral density induces a generally non- 
Mar kovian dissipative dynamics. 

As a main finding of this work, we have pointed out 
that the decoherence process happens in two stages. In 
the beginning, the purity decays like a Gaussian, while 
subsequently, Markovian decay sets in. If the qubit- 
oscillator coupling is strong or if the oscillator is strongly 
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driven, the major part of the coherence decays already 
during the first stage such that the relevant dynamics 
possesses a Gaussian time profile. Thus, with the trend 
towards ultra-str ong coupling between a qubit and a har- 
monic mode [2ll424| . Gaussian decay should become in- 
creasingly relevant. In the opposite limit of weak cou- 
pling, the first stage reduces the qubit coherence only 
by a small amount, rendering the relevant decoherence 
Markovian. The dephasing times in the two regimes ex- 
hibit distinct parameter dependencies, which we have 
determined analytically. Remarkably, in the Gaussian 
regime, the coherence time may be significantly longer 
than what one would expect from an extrapolation of 
the Markovian result. 

For a numerical description of the dynamics, we have 
solved the full Bloch-Rcdficld master equation for the 
qubit coupled to the driven oscillator in the frame of the 
original Hamiltonian. This has allowed us to obtain nu- 
merical results that are fully independent of the analyt- 
ical treatment. They have confirmed our predictions for 
the partially non-Markovian purity decay for the case of 
linear qubit-oscillator coupling. For quadratic coupling, 
the decoherence process is more involved. Nevertheless, 
it has been possible to obtain an analytical expression for 
the decoherence rate in the Markovian limit. Beyond this 
limit, our numerical solution indicates that decoherence 
is non-Markovian provided that the oscillator dissipation 
is very weak. 

Finally, we are convinced that our results on the de- 
coherence induced by a resonantly driven oscillator will 
support the design of future experiments with dispersive 
readout and its generalizations. 



Hamiltonian ([3]), such that 



Pose = - iUoWa, Pose 



*7 r r l i 7 i ( ^0 \ r r n 

— [x,\p,p osc \ + \ - -coth^^-^J[a;, [x,p osc jj. 

(Al) 



2. Average position 

Owing to the linearity of the quantum Langevin equa- 
tion for the dissipative harmonic oscillator, its position 
expectation value x obeys the classical equation of mo- 
tion 

I + 7! + u>lx = u) F(i). (A2) 

The response x(t) to the driving is most conveniently 
obtained by a time convolution with the Green's function 
Git), where 



= — o ■ 



(A3) 



and the inhomogeneity F(t). For F(t) = Acos(ilt — tfo), 
the solution reads 

x(t) = \G(Cl)\Acos(nt-<p-(j)o), (A4) 

where the phase shift <j> is the argument of the Green's 
function. The corresponding semiclassical state is a co- 
herent state with mean photon number n = ^A 2 \G(il)\ 2 . 



Position correlation function 
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Appendix A: The driven dissipative harmonic 
oscillator 

1. Markovian master equation 

In the absence of the qubit and for sufficiently weak 
oscillator-bath coupling, the dissipative dynamics of the 
oscillator is well described within Bloch-Redfield the- 
ory [3, [r| . Since additive driving of the form does 
not influence the dissipative terms |2{|, it is possible to 
use the master equation for the undriven dissipative har- 
monic oscillator and to simply add the time-dependent 



In the Markovian limit for the undriven dissipative 
harmonic oscillator implied in the master equation (|Aip 
in the absence of Hp, the equilibrium position auto- 
correlation function 

S SxSx (t) = ~{5[x(t),5x(0)}+) (A5) 

can be computed by help of the quantum regression the- 
orem [32j . It essentially states that S$ x Sx obeys the same 
equation of motion as the average position. Thus, 

Sx Sx 

+ uj SsxSx = 0. (A6) 

The initial value is 

*■«■!«- s"*"^)- (A7) 

while for symmetric ordering, its time derivative at t = 
vanishes. Thus, we can express the solution as Fourier 
integral of the Green's function (IA3|) which we evaluate 
via Cauchy's theorem. In the weak damping limit consid- 
ered herein, 7 < wq, we neglect the dissipation induced 
frequency shift and continue with the approximation 

S S x8x(t) = icoth (_^) e -7*/2 cosM . (A8) 
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